Introduction to Hubbard Model and Exact Diagonalization 
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Hubbard model is an important model in theory of strongly correlated electron systems. In this 
contribution we introduce this model along with numerically exact method of diagonalization of the 
model. 
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I. INTRODUCTION 

Hubbard modeli and its variants constitute an im- 
portant research topic in theoretical condensed matter 
physics, particularly in the context of strongly corre- 
lated electron systems. Most of the many-body tech- 
niques commonly used in condensed matter physics can 
be learnt in this context. Also there are some theoretical 
tools and concepts which apply to this model only. 

There are already some monographs^i^, which can be 
used by experts, along with some text books^i^i^ which 
can be consulted for further details of the various meth- 
ods. In this contribution our aim is to provide a smooth 
introduction to the model and exact diagonalization tech- 
nique used in dealing with Hubbard model. This work is 
based on the set of lectures given by the author in the first 
lUT schooU^ on strongly correlated electron systems. 

Analytical methods of solving the Hubbard model are 
all approximate, except in ID, where the so called Bethe 
ansatz provides an exact solution^. On the other hand 
there are exact numerical techniques, which are however, 
either computer time expensive or memory expensive. 
Therefore one is limited to rather small cluster sizes. 

A popular method to solve Hubbard model and many 
other models in condensed matter physics is exact diag- 
onalization (ED) of the models for small clusters which 
we will study at length in this set of lectures. We en- 
courage the reader to implement the method presented 
here in a simple f ortran program. In numerically exact 
diagonalization method one gets the 'exact' results at a 
high price, namely limitation to very small cluster sizes 
(about 18 sites for Hubbard model at half filling), which 
is essentially due to limited amount of computer RAM 
one can typically have. If one accepts some error bars in 
numerical results (which can however be systematically 
improved) , then a family of the so called Monte Carlo 
methods are methods of choice^. These methods are es- 
sentially exact. The accuracy of the results depends on 
how much computer time one would likes to spend. In 
this sense, these family of methods are time expensive, 
while ED is memory expensive. 

One of the important methods to deal with almost any 
model is Mean Field Theory (MET). MET ignores quan- 
tum fluctuations; hence becoming less accurate in lower 
spatial dimensions. Despite this, mean field treatment 
reveals the wealth of various condensed matter phases 
can emerge from a simple Hubbard model^. In MET one 



ignores both spatial and temporal fluctuations. It is pos- 
sible to retain the temporal fluctuations by performing 
a full fledged quantum dynamics of the problem. This 
is the subject of the so called 'dynamical' mean field 
theory (DMET>2.. 

There arc also a class of approximate analytic methods 
known as auxiliary particle or slave particle meth- 
ods, which are invented to deal with the large U limit 
of the Hubbard model™. This type of techniques are re- 
lated to the so called Gutzwiller projection which are 
devised to obtain approximate ground state of the Hub- 
bard model at half filling. Generalization of this method 
to deal with excited states is also there in the marketii. 

In the limit of large U, the charge fluctuations in the 
Hubbard model are frozen and only the spin of electrons 
can fluctuate. Thereby reducing the physics of the Hub- 
bard model to spin physics described by the so called t-J 
Modeli. 



II. HUBBARD MODEL 

A. Non-interacting electrons 

The Hamiltonian of a system of non-interacting 
fermions on a lattice of L sites labeled by i,j, etc can 
be represented in second quantization by 



Ho = ^Uj4cj 



where c'Acj) creates (annihilates) a fermion in a single- 
particle orbital (pj localized at site j. In condensed 
matter applications one can assume 0j's to be Wan- 
nier wave functions (Eouricr transform of Bloch orbitals). 
Eermionic operators satisfy the anti- commutation rela- 
tions, 



{cl,Cj} — Sij, others = 0. 



(1) 



The coefficients tij characterize the single-particle matrix 
elements 



2y72 
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V \(j)j 



dxMx) [ -^^ + ^'(^) ) '/'j(^) (2) 



For many practical purposes it suffices to assume that tij 
is none-zero, only when i,j are nearest neighbors in which 
case it is usually denoted by —t, so that the Hamiltonian 
written in manifestly hermitian format becomes 






c]c. 



(3) 



Assuming the periodic boundary conditions (PBC), the 
Hamiltonian of the system will be invariant under trans- 
lation. The irreducible representations of the translation 
group (due to abelian structure of the group), are one- 
dimensional (i.e. numbers of type e'^). Hence the one- 
particle (or non-interacting=free) Hamiltonian ^ can be 
diagonalized by a Fourier transformation [9 <-^ kj) 
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Then the Hamiltonian in 



basis becomes 



^o = E' 



:fcc[cfe, efe = -2icosfc. 



(4) 



(5) 



where eu determines a cosine dispersion and represents a 
band structure with L allowed k values in the first Bril- 
louin zone. Of course a simple cosine may not be a good 
approximation of the realistic band structure of solids 
(spaghetti plots). To mimic the realistic band structures, 
one can add further neighbors' hoppings which generate 
higher harmonics of the simple cosine band. 

ID tight-binding band dispersion 




FIG. 1: A tight-binding band dispersion in ID with electron- 
like Fermi surface. 



If the number N of the electrons is equal to number 
L of the sites, then each allowed k state can be occupied 
by two t and | spins. Hence the ground state of Hq is 
constructed by filling the lower half (e^ < 0) of the band 
dispersion of Fig. [T] which is denoted by circles on the fig- 
ure. Since half of the band is filled, the N = L situation 
is called half-filling^. This state is known as a Fermi sea 



state, usually denoted by |FS). 
sentation of this state is: 



Second quantized repre- 



|FS) 



n ^TclilO), 



(6) 



k<kp 



where |0) is vacuum state (empty lattice) and kp is the 
largest occupied k value (7r/2 here), known as Fermi 
wave-vector. This expression in first quantized notation 
corresponds to an slater determinant. So the ground 
state of a non-interacting Hamiltonian Hq is character- 
ized by a single slater determinant. 

This state is an eigen state of the Hamiltonian 



Ho|FS)=£;o|FS), 
where Eq is the total energy given by 

'•^/^ dk 



(7) 



Ef ' dk 
eknka= / 7^2 X (-2icosfc), (8) 

111^ /9 J-tt/2 27r 

|/c|<7r/2,a- ' 

where the single-particle energy summation is carried out 
over the occupied k states of Fig. [T] shown by circles. 
Obviously the depicted state has the lowest energy, since 
the k states have been filled such that lower energy single- 
particle states are filled first. In this spirit and excited 
state is create as follows: leave on state k with |fc| < 7r/2 
empty, promoting its electron to another state k-\-q such 
that |fc + (j| > 7r/2. This excited with a hole left behind 
in state k and an electron created in state k + q with, say 
t spin is called a 'particle-hole' excitation: 

TJolV-n) = i/o4.+,TCfeT|FS) - <,|^„) (9) 

which is again an eigen state of Hq with energy e^ {q) = 
Sk+q — Sk- This excitation carries a center of mass mo- 
mentum q. The above discussion can be straightfor- 
wardly generalized to higher dimensions. 

Evidently the half-filled state |FS) characterizes a 
metal, as the energy of the particle-hole excitation can be 
made arbitrarily small. The above band picture always 
gives a metal for an odd (in this case one) number of 
electrons per unit cell. However, as we will argue in the 
following, there might be situations in which this simple 
prediction of the band theory fails drastically. 



B. Electron-electron interaction 

When we have only Hq (Eq. ^ in the Hamiltonian, 
then the minimization of energy is achieved by filling the 
k states independently with electrons of opposite spins. 
Such k space picture in real space translates to equal 
probabilities p = 1/4 for four possible occupation of a 
single site: 



o 



(10) 



The most general form of interaction in second quanti- 
zation representation which can be added to Hq is of the 
form 






(11) 



where a = {a, i} is a collective name for the site index 
I, and spin index a the two-particle matrix elements is 
given by 

V/3a = f dxdx'rjx)r,ix')Vi\x - X'\)MX)M^') 



(12) 
usually in metals with appreciable density of states D{ep) 
at the Fermi level, the Coulomb potential y(|x — x'| — r) 
is screened and obtains the form 



^ — rk'T 



Vir) 



(13) 



This limit known as band limit always describes a metal. 
The U term on the other hand increases the cost of charge 
fluctuation, leading to an insulator in the opposite limit 
[/ — > oo. Therefore there must be a critical value Uc of 
the order of band width W = 2zt (where z is the coordi- 
nation number), beyond which one has an insulator. This 
(first oder at zero temperature) phase transition is known 
as Mott metal-insulator transition (MIT). The only tech- 
nique which can handle this model at arbitrary values of 
U and for arbitrary filling is ED, which will be described 
in next section. 



III. EXACT DIAGONALIZATION 

The easiest way to describe the essence of this method 
is by the example of a two-site Hubbard model. This 
toy model consists in two sites labeled 0, 1. In this case 
the Hubbard model written explicitly (in units in which 
t = 1) reads: 



H 



= Ht + Hu 



where k^^ is the so called screening length. For the d 
electron systems where the overlap between the atomic 
wave functions is small (smaller t) , one has narrower band 
which is synonymous to larger DOS at the Fermi level 
the screening length is usually on the scale of the Bohr 
radius a^. Therefore the most important term among 
all possible iJLv(3a matrix elements is when all the indices 
correspond to the same site j. In such case, the Pauli 
principle forces fi = a =t and j3 = v =J,. Denoting the 
corresponding matrix element with — 2C/ one gets for the 
screened interaction 

y — — [/\^c^ c^ c tc. I = U'S^ nin-\ . (14) where the site index increases from right to left. 



(16) 



where sites are labeled as follows: 



If we add this term to Hq, we obtain the celebrated 
Hubbard model 



H 






^ia^j'^ 



Cj^C^. 



U^n^^riji. (15) 



For the Hamiltonian in [/ — > oo limit the doubly occu- 
pied configuration of a single site in Eq. pop is going to 
cost a large energy Utij^tiji = t/(l)(l) for each doubly 
occupied site. Therefore presence of such term violates 
the equi- probability of four possible states of a single-site, 
thereby inducing some kind of correlation by minimizing 
the double occupancy. Therefore at half-filling, and for 
large U the ground state charge distribution adjust itself 
to avoid doubly occupied sites as much as possible; i.e. 
each site is occupied by a single electron, and that more- 
over, large value of U makes charge fiuctuations around 
rij^ = 1 configuration very expensive. Therefore charge 
fluctuations are frozen and one has an insulator known 
as Mott insulator. 

For finite values of U, the two terms of the Hamil- 
tonian (jlSp compete with each other. The kinetic en- 
ergy term (corresponding to f/ = 0) tends to delocalize 
electrons by putting individual electrons in Bloch states. 



A. Organizing the Hilbert space 

To organize the Fock space for this Hamiltonian, one 
first notes that the number operator N = '^,j^ rij^ com- 
mutes with the Hamiltonian. Therefore one can consider 
only Hilbert space corresponding to a fixed value of TV. 
Let us consider N ~ 2 for this toy model which corre- 
sponds to half-filling condition. The next question arises 
with regard to the total spin of the electrons: whether 
they are 11, H or tj? First two cases represent triplet 
state, while the last one corresponds to 5" = (more pre- 
cisely {11 — 11) is a singlet). Formally it can be checked 
that the total Sz ~ 1/2 X); ('^jT ~ ^ji) ^l^*-* commutes 
with the Hamiltonian and hence a conserved quantity. 
Therefore there would be no matrix element of the Hamil- 
tonian connecting sections of the Hamiltonian with dif- 
ferent values of Sz- The structure of the Hamiltonian 
will be block diagonal where each block corresponds to a 
fixed value of Sz- To see this block-diagonal structure, we 
confine ourselves to TV = 2 with both triplet and singlet 
spins. 

In sector with quantum numbers N = 2 and Sz — 0, 
Hilbert space is six dimensional with six possible basis 



states |(/)j) with J = 1 ... 6. 





beginner 


expert 


computer 


J 


algebraic 


picture 


binary /| J| / 


1 




T T 


|00)J11)| 3 3 


2 


t t IX 

'^ox'^otI/ 


OTi 


|01)i|01)T 1 1 5 


3 




T i 


|01)x|10)t 1 2 6 


4 




i T 


|10)x|01)t 2 1 9 


5 


t t i\ 


TiO 


|10)x|10)t 2 2 10 


6 


t 1 1\ 


i i 


|ll)x|00)| 3 12 



(17) 



First row indicates three major sets of cohimns: 

• The second column labeled 'beginner' is a way a 
beginner would define and work with a set of 6 basis 
states \(j)i} to 106)- 

• The third column labeled 'expert' is a way an ex- 
pert works with these basis states. 

• Columns number 4 to 7 are related to the way a 
computer organizes and works with the basis states. 

Explanation of each columns is as follows: 

1. In the first column there is an integer J = 1 ... 6 
which labels the basis states in the Hilbert space. 

2. Second column shows how to obtain the configura- 
tion depicted in third column by acting with cre- 
ation operators on the vacuum |). As an example 
look at the first basis \(j)i): Physically, it describes 
two t electrons in sites number 0, 1. Therefore we 
have two choices: 



10?) 



"it"^ot I/' 



''Ot''1T 



These two choices by Fermionic anti commutation 
relations are negative of each other: 10°) = — |0i). 
Then the question is: which one is correct way of 
representing |0i)? The answer is that, both of them 
are fine, the same way one can describe ordinary 3 
dimensional vector space with basis ei = x, 82 = 
y, 63 = z, or say, ei = -i, 62 = y, 83 = z. The 
important point is to stick to one convention during 
the entire matrix and vector manipulations. For 
example, to construct the above table we choose 
the following convention: (i) | spin operators sit to 
the right of | spins, (ii) The order of site indices 
increases from right to left. 

3. Third colunni is a pictorial representation of the 
basis vector in real space. We will learn in the 
following how to work with this intuitive represen- 
tation of the basis states. 

4. The fourth column represents the states as direct 
product of spin-J, with spin-]" state states, where 
the occupations in up and down spin sectors are 



represented by sequence of 0, 1 bits giving rise to a 
binary representation of the basis states. 

For models like, Hubbard model where t and | 
spins do not admix, this separation makes a good 
sense. To see this, consider the effect of a term like 

J 



chcjT^ on a typical state 



<i- 


■^iNii tS'lT ■• -S'mT 


<.r 


■4ji)i[4T---4 



jmtjDt 



which can be written as 

= \c^. , ...c^ ,] (-l)^(-l)^c^,c,t[c^ , 



h4i 



"■iiviJL-jit 



"iMTJ 



^jmTJ 



(18) 



Here two (—1)^ factors arise from moving each of 
the cl operators through a length iV sequence of c! 
operators. Therefore, when operating in spin-\ sec- 
tor, we need not worry about spin-l configuration 
and vice versa. 

5. The last three column indicate the way a computer 
understands and stores these basis state. There is 
only one integer / stored on computer which fully 
represents the occupation pattern of the "f and | 
electrons when transformed to binary representa- 
tion of fourth column in Eq. p7p . Given /, one can 
obviously find out /| and I^ through the relation 



/ = 2^/x + /t 



(19) 



and vice versa. First and last columns of ([T7)) do 
actually tabulates an array / = T{J) which for any 
given J, returns the corresponding I. The value of 
/ fully specifies the state. One therefore only to 
extract the bits of / and find a way to work with 
bits of integer /. 

There remains only a final note on how we have labeled 
the states. We have labeled the configuration represented 
by / = 9 as fourth (J = 4) basis vector, etc. We could 
have labeled them in any order, so that 7 = 9 would 
have corresponded to, e.g. first (J = 1) basis vector. 
The above convention which was suggested by Lin and 
coworkersi^, has the advantage that the table T can be 
searched given a value for / in a fast way in order to 
find out corresponding J value. The essential idea of 
this convention is the foUowingi^: For any sector you are 
interested in, just choose the labels J in such a way that 
when the above table T is constructed, the / values are 
ascending function of J. 



B. Acting with operators on the basis states 

If the action of an operator on a complete basis set 
is known, then the operator is completely specified. 



Eq. ([T6|) has two type of terms: Hu and Ht . When H{j 
acts on a basis states, it gives non zero contribution for 
each site j in which both riji and riji are 1. Therefore 
the effect of Hjj on any basis state gives the same state 
multiphed by the number of doubly occupied sites xU. 
Hence the effect of Hu on {(pi), i = 1,3,4,6 is zero and 
on \(f)2) and j^s) is just U. The Hu part is diagonal in 
occupation number representation: 



Hjj^U 



/O 0\ 

10 





10 

Vo 0/ 



(20) 



Now let us concentrate on the effect of Ht on our basis 
states. Consider, e.g. the sate |02). Ht in Eq. (fT6|) has 
4 terms. But only 2 of them give non zero contribution 
when they act on 102): The one which allows an ] spin to 
hop from site to 1, i.e. cLco|, and another one which 
allows the [ spin at site to jump to site 1, i.e. cji cq^. 



Ht\cl>o) 



-t {c\^co^ + c\^caij cl^cl^\) 



= -t 



-1| Cot 






Cl^CO^CgJ^ Cq| 



t J 



t f-cT^cJ^coTcJ-rl) +c{^cJ^coicJ^ 




-t 



Ci|CqJ) + C;^|Cg|| 



where in the third line we have used the Fermionic 
anti-commutation relations with associated minus signs 
needed for each exchange of Fermionic operators with dif- 
ferent indices. Also in the last step we have used the fact 
that Cjo-c]^|) = (1 — nja)\) = |), as the vacuum state |) 
contains no particles. Rearranging the Fermionic opera- 
tors to comply with our convention we get 



site 1, or to move a | spin to site 1. In the former case 
one gets (— l)°|(/i3) where the exponent is because the t 
spin at site passes through no other t spin when it hops 
to site 1. Similarly the later case gives [(^4). According 
to the Hamiltonian ([16]) , each of these processes happens 
with an amplitude —t and thus one can pictorially see 
that 

i?t|</'2)=-i( + |03) + |04)). 

With this in mind it is almost trivial to check that 



i?d01> = 


= 


i/t|02) = 


= -t(|03) + |04)) 


i?*|03) = 


---t{\h) + \cp5)) 


Ht\(j>A) = 


= -t(|<^2) + |05)) 


i?t|05) = 


.-t (103) + 104)) 


i?t|06) = 


= 



(22) 



Therefore in this basis the hopping term has the following 
matrix representation: 



/O 0\ 
110 
10 10 
10 10 
110 

Vo 0/ 



Ht^-t 



(23) 



Therefore the matrix representation of the entire Hamil- 
tonian in this basis becomes 



H = 



/o 





0\ 





U -t -t Q 








-to -t 








-to -t 








-t -t u 





\o 





0/ 



(24) 



ift|02) = -t(+|03) + |04)) 



(21) 



Clearly the above block diagonal structure has a parallel 
with the following form for the Sz matrix: 



This way of working with commutations and algebra is 
not a convenient for putting on computers. As we already 
showed in Eq. (fT5)) . the up and down spins hop separately. 
Therefore, when an spin "f electron hops from site say j 
to site i, one only needs to count a (—1) factor for each 
"f spin electron over which it passes. Since when Cj| has 
to moved through a chain of c-,-^ operators with j' ^ j. 



until it reaches cL, there they form a (1 — n^i) operator 

which commutes with all other remaining c-,,-> operators 
to reach the vacuum |) where it produces |) itself. 

With this argument in mind, one can most conve- 
niently work with the pictorial representation of states: 
The effect of Ht on |02) is to either move an | spin to 



S, 



/+1 




\ 













I 




-1/ 



(25) 



In other words state |0i) has a quantum number S^ = 
-1-1, while 106) has S^. = —1. The set of states |02) . . . |05) 
belong to Sz = sector. If we had confined ourselves to 
Sz = sector, we would have obtained a four dimen- 
sional Hilbert space in which the Hamiltonian could be 
represented by the 4x4 block of Eq. ([M|) . 



1 1 from a 



The method deseribed after Eq. pT|) has the advantage 
that for large matriees, precisely the same steps can be 
taken with a simple computer code that implements the 
logic described above. When coding the procedure we 
have the binary representation in column 4 of Eq. (|17p 
in mind for the occupation pattern of ] and | electrons. 
But there is only one integer / (7th column) stored on 
computer the binary representation of which precisely 
corresponds to column 4. Extracting /| and 
given / in f ortran is as simple as 

I_up=mod(I,2**L) ; I_dn=I/2**L 

Most of the standard high level programming languages 
such as C++ and f ortran 90 have appropriate intrinsic 
functions for bitwise operations on integer which are able 
to access and examine or change individual bits of a given 
integer p (can be /| or /^). Therefore an integer I com- 
pletely specifies the occupation pattern in each spin sec- 
tor. For example the Intel Fortran compiler has following 
commands for bitwise manipulation of an integer p 



summarized as follows 



IBSET(p,b) 
IBCLR(p,b) 
BTEST(p,b) 



(26) 



For more details please consult the Intel Fortran language 
manuali^. Also more implementation notes can be found 
in Lin et.alJ^. 



C. Diagonalization of the matrix: Space 
symmetries 

So far we have made use of the symmetries of the 
Hamiltonian itself to reduce the dimension of the ma- 
trix to be diagonalized. For example conservation of Sz 
(equivalent to [Sz,H] = 0) reduces the 6x6 matrix to 
two 1x1 (trivial) and one 4x4 matrix. The commuta- 
tion relation between e.g. Sz and H is quite general, and 
independent of the geometry of the lattice in use. 

Now we would like to make use of the spatial symme- 
tries in order to diagonalize the 4x4 matrix in Sz ^ 
sector manually. For a two site problem composed of 
sites 0, 1, there is a mirror reflection operator M which 
has the following action on real lattice: 



A/(0) = 1, M(l) = 



(27) 



Corresponding to operation M which is a member of spa- 
tial symmetry group of the underlying lattice, there is an 
operator A4 in Hilbert space which acts on the state vec- 
tors. For example consider the effect of A4 on state |03) 
which reads 

Mclcl^lO) = 4/(o)i4/(i)TlO) = ^i^tlO)' (28) 
or A^|03) — \(f)4). Similarly, A^|04) — j^s), etc. which is 



X|0l)=- 


-\4>i 


M\(j)2) = 


\4>5 


M|03) = 


104 


A^|</'4) = 


103 


M\^5)^ 


102 


Mlcj^e)^- 


-|06 



The first and last lines of the above equation indicate that 
states \<f)j) with j = 1, 6 are eigen states of Al (parity op- 
erator) with eigen values ±1. Line numbers 3,4 however, 
indicate that states |03) and |04) do not have definite par- 
ity. Instead a new combination |0±) = -i=(|03) ± |04)), 

has a definite parity: M\(f>±) = ±|0±). Similarly from 
102) and 105), symmetric and antisymmetric combina- 
tions will have definite parity. 

The eigenvalues ± of parity operator M partition the 
Hilbert space into two pieces which belong either to +1 
eigen value: 



(29) 



M {\(l>3) + \cj)4)) = 
A^ (102) + 105)) = 


= + (103) + 104)) 
= + (102) + 105)), 


en value: 




A^|0i) = -|0i) 

A^|06) = -|06) 
M (103) - 104)) = 
M (102) - 105)) = 


= -(103) -104)) 
= -(102) -105)) 



(30) 



If we knew this when finding out the effect of Ht in 
Eq. p2)) . we would have acted on the following states: 



l^i) 

1^3) 

1^4) 

IV'5) 

IV'e) 



\4'i) 

(102) 
(103) 
(103) 
(102) 



14',)) 
104)) 
104)) 
105)) 



(31) 



Eq. (PT|) defines the so called, symmetry adopted basis in 
which the action of, say, Ht is enormously simplified: 



i?*|V'i) = 


= 




Ht\i^2) -- 


= - 


2<|^3 


Htm = 


= - 


2i|V'2 


HtliJi) = 


= 




Htm = 


= 




Ht\A.) - 


= 





(32) 



Similarly the only none zero matrix elements of Hu are 
given by: 



Hum=U\ij2), Hum^Um 



(33) 



How do we generate a symmetry adopted basis? There 
is a very powerful theorem in group representation theory 
which in case of one dimensional representations is easy 
to implement and reads^^: 



^ 



(p) 



^rp(i?)*i?0 



(34) 



In this equation (j) is an arbitrary state to begin with. R is 
an element of the symmetry group which in this case and 
be either X or A^. Fp is the irreducible representation 
of the group which for group {X, Al} are numbers ±1 
for even and odd parity. Reader can easily check that 
feeding = \4>i), with i = 1, . . . , 6 is going to generate 
the symmetry adopted basis Eq. (pij) . Similar technique 
can be used to generate states with definite k values when 
dealing with problems of translational invariance. 

In symmetry adopted basis total Hamiltonian in Sz ~ 
±1 sectors ({|'!/'i), IV'e)} sub-space) remains diagonal with 
eigen values equal to 0, as it was. In 5*2 = sector it will 
become: 



H = 



I ^ 


-2i 


\ 


-2t 

















V 





u J 



(35) 



In this sector as well, there are one eigen value equal to 
which corresponding to eigen state |V'4)- Note that \i(;4) is 
the Sz ~ component of a triplet. Therefore the triplet 
state lies at zero energy: 



Et = 0, for triplet state (5*^ = +1,0, -1) 



(36) 



The fact that {|0i), l^s), |<^6)} form a triplet is consistent 
with having odd spatial parity in Eq. (pO|) . 

The eigen value corresponding to {tp^) is U which is 
always positive and above the triplet Et ~ 0. 

To find out the remaining two eigen values in the 
{IV'2), I "03)} sector, we note that first of all, this sector 
has even spatial parity, and hence is spin singlet. Hamil- 
tonian is U/22 + 11/2(7 2 ^ 2icra;, where (Jx,o'z a-re Pauli 



matrices, and T is the unit 2x2 matrix, so that the 
eigenvalues at spin singlet sector become 



Ef = U/2 ± ^/{U/2)^ + 4i2 



(37) 



Therefore the ground state is a singlet with energy E~ , 
while the first excited state is a triplet with Et = 0. 
Second excited state is |'05), with energy U, and final 
excited state has energy E^ . 



IV. STRONG CORRELATIONS AND SPIN 
PHYSICS 

As we already saw in previous section, the ground state 
of the two site Hubbard m odel is a singlet with energy 
E^ = U/2 - y^U'^/4: + W^, The ground state wave func- 
tion 



\E7) = 



) = 4|V'2) + (u + Vu^ + le) IV3) 



(38) 



in large U limit is dominated by lip^) ~ j^a) + |04) in 
which there is no doubly occupied configuration, and 
hence charge fluctuations are suppressed. 

Since the first excited state is at Et = 0. The sphtting 
between these two states for large U ^ t is: 



J = Es-Et^ U/2 - ^JlPjI+W 



4t2 



(39) 



Therefore the singlet state is slightly below {—At^/U) the 
triplet state. This indicates that in large U limit, the low- 
energy physics of Hubbard model is given by spin fluc- 
tuations which are anti ferromagnetic (singlet has lower 
energy). This observation in a two site Hubbard model is 
indeed very general and it can be shown using a unitary 
transformation that the Hubbard model at large U limit 
can be mapped into the so called t-J model, where there 
are AF spin fluctuations along with hoppings restricted 
to subspace with no double occupancy^. 
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